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ABSTRACT 

The  interaction  of  three  jets  is  studied  both  theoretically  and 
experimentally o   The  angle  of  deflection  and  the  velocity  and  turbulence 
distributions  were  determined  through  the  use  of  a  hot-wire  anemometer. 
The  free -streamline  theory  has  been  applied  to  the  determination  of  the 
jet  deflection  and  the  results  so  obtained  were  found  to  be  in  good 
agreement  with  those  obtained  experimentally „   The  study  is  a  convincing 
evidence  of  the  fact  that  a  proper  combination  of  inviscid  flow  analysis 
together  with  experimental  results  can  provide  the  necessary  guidance 
for  the  design  of  proportional  amplifiers. 
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I.   INTRODUCTION 

The  interaction  of  compressible  or  incompressible  jets  has  been 
subject  to  study  for  at  least  the  past  fifty  years  for  the  purpose  of 
predicting  the  velocity  and  turbulence  distributions  in  the  combined 
jets  and  in  general  for  the  purpose  of  predicting  the  performance  char- 
acteristics of  all  sorts  of  machinery  where  such  jet  interaction  plays 
a  significant  role.   Most  of  the  theoretical  analyses  have  been  restricted 
to  the  study  of  the  interaction  of  jets  of  an  inviscid  fluid,  and  the 
complex  effects  of  viscosity  and  the  effects  of  the  bounding  walls  have 
been  largely  ignored. 

The  advance  of  the  fluidics  technology  and  the  discovery,  among 
others,  of  proportional  amplifiers  have  renewed  a  strong  interest  in 
both  the  theoretical  and  experimental  investigation  of  jet  interaction. 
Proportional  amplifiers  were  designed  and  operated  either  through  the 
use  of  the  "pressure  control"  concept  or  through  the  use  of  "momentum 
control"  concept.   In  either  case,  both  the  single-sided  and  double- 
sided  devices  were  considered.   The  pressure  controlled  device  is 
operated  on  the  basis  of  the  differential  pressure  between  the  two 
cavities  placed  on  the  two  sides  of  the  resultant  jet.   The  vents  are 
not  open  to  the  atmosphere.   The  momentum  controlled  proportional  ampli- 
fier operates,  as  its  name  implies,  on  the  momentum  balance  of  the  com- 
bining jets.   Thus,  it  is  apparent  that,  whereas  the  pressure-controlled 
amplifier  depends  primarily  on  pressure  to  deflect  the  jet,  the  momentum- 
controlled  amplifier  depends  primarily  on  momentum.   For  certain  con- 
figurations where  the  venting  of  the  side  cavities  is  not  complete  or 
where  the  entrainment  from  the  two  sides  of  the  jet  is  not  completely 
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attainable,  then  it  is  necessary  to  take  both  pressure  and  momentum  ef- 
fects into  account  in  attempting  to  design  amplifiers  optimized  for 
particular  conditions.   It  should,  in  passing,  be  noted  that  the  use 
of  both  control  mechanisms  tends  to  produce  devices  whose  band-widths 
are  relatively  small,  even  though  they  may  have  significantly  higher 
pressure  gains.   We  shall  confine  our  attention  here  to  a  thorough  study 
of  the  momentum-controlled  devices  or  beam-deflection  amplifiers,  both 
for  the  purpose  of  understanding  their  fluid  mechanics  and  for  the  pur- 
pose of  contributing  to  a  predictable  design  of  proportional  amplifiers. 

The  literature  concerning  the  fundamental  behavior  of  beam-deflec- 
tion devices  is  still  relatively  limited  but  is  growing  rapidly  with  the 
increasing  realization  of  the  design  limitations  of  current  fluidics 
technology o 

The  prediction  of  the  performance  of  such  devices  depends  on  a 
knowledge  of  the  exact  mode  of  interaction  of  the  jets.   There  are, 
obviously,  many  possible  geometric  and  kinematic  variables  for  the 
interaction  process  and  it  is  not  economically  feasible  to  experimental- 
ly investigate  the  characteristics  of  the  mode  of  operation  of  each  and 
every  possible  combination.   This  realization  necessarily  led,  in  the 
first  decade  of  fluidics,  to  the  evolution  of  several  flow  models  with 
varying  degrees  of  approximations  and  complexities.   There  are,  to  be 
sure,  essential  differences,  as  well  as  similarities9  among  these  models, 
and  all  must  be  carefully  assessed  before  a  design  choice  is  made. 

Dexter  [Ref„  1]  employed  the  "submerged  jet"  model,  assuming  that 
viscous  mixing  with  the  surrounding  fluid  has  reached  the  stage  where 
the  velocity  profile  of  the  deflected  jet  is  of  Gaussian  form.   Moynihan 
and  Reilly  [Ref.  2]  employed  the  "self  preservation"  model  in  which  each 
jet  is  assumed  to  have  retained  its  identity  during  interaction  and  is 

14 


not  subjected  to  any  viscous  effects  as  it  proceeds  downstream,   Douglas 
and  Neve  [Ref.  3]  extended  the  works  in  the  references  just  cited  by 
applying  the  two  models  to  the  appropriate  regions  of  the  resultant 
jet. 

The  present  work  presents  the  theoretical  results  obtained  through 
the  use  of  the  free-streamline  theory  for  symmetric  control-jet  con- 
figurations and  the  experimental  results  obtained  with  a  hot-wire 
anemometer. 
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II.   THEORETICAL  CONSIDERATIONS 

As  in  all  applications  of  the  free-streamline  theory,  the  pressure 
and  magnitude  of  velocity  along  the  free-stream  lines  are  assumed  con- 
stant.  It  is  further  assumed  that  both  the  pressure  and  velocity  are 
continuous  along  the  separation  lines  between  the  power  jet  and  the 
control  jets. 

Experiments  show  that  each  control  jet  retains  a  large  measure  of 
its  integrity,  and  the  power  jet  undergoes  an  acceleration  or  a  "vena 
contracta"  effect  during  the  initial  face  of  the  interaction.   The 
combined  jet  attains  a  Gaussian  profile  only  after  a  distance  of  two 
to  six  nozzle  widths  from  the  pivoting  point,  depending  on  the  flow 
ratio  and  the  geometry  of  the  system. 

It  is  apparent  from  the  foregoing  that  the  flow  in  the  region  in 
which  the  momentum-exchange  and  the  turbulent  mixing  take  place  does 
not,  strictly  speaking,  follow  the  assumptions  made  in  the  free-stream- 
line theory.   The  complicated  interaction  between  the  jets,  and  the 
pressure  of  turbulent  mixing,  make  the  analytical  predictions  of 
three-dimensional  flow  field  extremely  difficult.   To  obtain  analytical 
solutions,  a  simplified  incompressible  flow  model,  which  retains  some 
of  the  major  features  of  the  actual  flow  field,  must  be  employed.   Al- 
though the  probable  error  introduced  through  the  use  of  such  a  model  is 
very  difficult  to  estimate,  previous  analyses  place  a  reasonable  degree 
of  confidence  in  the  applicability  of  the  model  to  certain  features  of 
the  jet  interaction  phenomenon. 

A.    TRANSFORMATION  PLANES 

A  schematic  drawing  of  the  flow  configuration,  together  with  the 
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resultant  jet,  is  shown  in  Fig.  la.   This  configuration  may  be  trans- 
formed to  an  iTL-plane  through  the  use  of  the  Planck's  transformation 
with  the  usual  assumptions  of  the  free~streamline  theory.   Thus,  writing 

H„   <tv\^   *  ce  n l) 

and  carrying  out  the  necessary  manipulations,  one  obtains  Fig.  lb. 

It  is  a  well-known  fact  that  any  polygon,  such  as  the  one  shown  in 
Fig.  lb,  may  be  transformed  to  the  upper  or  the  lower  half  of  a  real 
plane  through  the  use  of  Schwartz-Christof fel  transformation  given  by 

where  M  and  N  are   constants   to  be   determined   and    Zi  and    ^i   are   given 
by  £,  =  o  rtal  £*=■-<         r4^k         £r  =  -^ 

Then  one  has 


r 


n  = 


Mtdfc  uU  (3) 


+  Kl 


The  evaluation  of  this  integral  is  shown  in  detail  in  Appendix  I. 


The  result  is  that 


-k 


(4) 


(5) 


(6) 
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Table  I  gives  the  corresponding  values  between  the  points  in  the  -Q_  -plane 
and  the  t-plane  shown  in  Fig.  lc. 

The  relationships  between  the  parameters  shown  in  the  t-plane  and 
various  velocity  ratios  are  obtained  as  follows: 

Through  the  use  of  the  conditions  at  point  A 


one  has 


_  VA  W 


n  _  VM  _       *> 

v.  vu      \.+  /TIF 

Through  the  use  of  the  conditions  at  point  C 

one  has  ^J  ^    \/c*-kz 


Through  the  use  of  the  conditions  at  point  G 

one  has  VJ  V  -  ^^-k7- 


t 


Through  the  use  of  the  conditions  at  point  E 


one  has 


Sin  A  *  .1  f2-!  Qnc|  COi  ft  s  >i    '-k^ 

1  i2-e  I       Vf\^ 


4-^ 


(7) 


(8) 


(9) 


(10) 
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Finally,  through  the  use  of  the  conditions  on  the  free  streamline,  we 
obtain  /    . 


or 


J2.  ^  L  B 


±& 


Si. 

e   =, 


\ 


-  k 


tx-k' 


<• 


e-k2- 


COS.  ©  +  U  Surv  © 


(11) 


It  should  be  noted  that  the  conservation  of  mass  yields 


d  =  (of+#)b+ooL 


The  t-plane  is  comprised  of  three  sources  and  one  sink.   Thus,  the 
complex  potential  for  the  flow  in  the  upper  half  of  the  t-plane  may  be 
written  in  terms  of  the  known  potentials  of  sources  and  sinks  as: 

w=  vj'i  Wt^Vv*^  ^fc  -VC!s  Wt+cN-  V€»k  ^[t-l)  (12) 

"C       T  '    X         X   v    '  it     °' 

where  W  -  ty  +  ^  j^     and     <2>     is  the  complex  potential  and  j 
the  stream  function. 

Taking  the  derivative  of  w  with  respect  to  t,  making  use  of  Eqs. 

(7),  (8),  and  (9),  and  the  conservation  of  mass  relationship,  one  has 


dfc 


:  ^ 


_  -^_   +  1  £•  +    ^  ^    -v  iLk 


^  L    t+4 


t-vc       t- 


JJ 


(13) 


The  distances  between  various  points  in  the  z-plane,  the  contraction 
coefficient,  and  the  angle  of  deflection  may  be  evaluated  as  follows: 


5?s  ^[jd*] 


(14) 
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and  noting  that  DF  =  a  +  2s9  one  has 

dw  /~\ 

Replacing  —  from  Eq.  (13)  and  SL  from  Eqs.  (4)  and  (6)  9   and  taking 

the  Cauchy  principal  value  of  the  integral  to  account  for  the  discon- 


tinuity at  t  ■  -f,  we  finally  have9 


+■  2.5  =  d  6i  n  fo  +  Uoo 


£-0 


r        f  i 


lunrv 
£-»0 
Q  -*oo 


Q-*oo 


J~ 
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The  details  of  the  integration  are  shown  in  Appendix  II.   The  results 
may  be  expressed  ass 
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Equation  (16)  establishes  the  relationship  between  the  various  jet 
velocities  and  the  amplifier  setback.   Another  expression  is  needed  to 
establish  the  symmetry  of  the  geometrical  configuration.   This  may  be 
done  by  performing  an  integration  from  D  to  F  along  the  free -streamline 
and  writing  the  imaginary  part  equal  to  zero  so  as  to  ensure  that  the 
two  sides  of  the  control  jets  are  along  the  same  line  and  that  there  is 
no  unsymmetry.   Thus,  writing 
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or 
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and  replacing  -r~     from  Eq.  (13)  and  ±L   from  Eqs»  (4)  and  (6),  one  has 


-I 


(is; 


Details  of  the  integration  are  shown  in  Appendix  III.   The  analysis 
yields: 
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where  l\C^)  is  the  complete  elliptic  integral  of  the  first  kind  and 

I  I?  i  °^Z >  ^  i    *s  tne  comPlete  elliptic  integral  of  the  third  kind- 

2 
Since  the  parameter^  of  the  elliptic  integral  of  the  third  kind 

2 
shown  in  Eq.  (19)  was  not  always  in  the  range  of  0^.  <5(  ^  1,  a  sub- 
routine program  was  written  to  reduce  all  the  complete  elliptic  inte- 
grals of  the  third  kind  into  the  sum  of  complete  and  incomplete  elliptic 
integrals  of  the  first  and  second  kind  for  which  there  are  tables  and 
subroutines  available  at  the  computer  facility  of  the  Naval  Postgraduate 
School. 

The  reductions  were  made  through  the  use  of  the  transformations 
given  in  Ref.  [4],  and  through  the  use  of  Heuman's  Lambda  and  Jacobi's 
Zeta  functions,  as  given  in  Ref.  [5].   The  aforementioned  reductions  are 
shown  in  Appendix  IV.   The  evaluations  of  the  resulting  expressions  for 
the  discontinuity  at  t  ■  -f  is  shown  in  Appendix  V.   The  results  obtained 
for  Jj.  ^  \  and  b  ^  |   through  the  use  of  the  expressions  given  above 

will  be  presented  together  with  the  experimental  data  in  connection  with 
the  discussion  of  results. 
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The  experimental  apparatus  consisted  of  a  jet  assembly,,  hot-wire 
aneomometer  system,  and  a  velocity  calibrator, 

A.    JET  ASSEMBLY 

As  shown  in  Fig.  2  and  3,  the  entire  assembly  consisted  of  three 
nozzles,  one  for  the  power  jet  and  two  for  the  control  jets.   Each  noz- 
zle was  1/4  inch  wide  and  2  inches  high.   This  resulted  in  an  aspect 
ratio  of  eight.   The  jet  assembly  was  sandwiched  between  two  1/4  inch 
plexiglass  plates. 

The  width  of  the  control  nozzle  was  adjustable  from  1/4  inch  to 
1/2  inch.   For  all  the  tests  reported  here,  the  control  and  power  jet 
nozzles,  as  well  as  the  setback,  were  maintained  at  widths  of  1/4  inch. 

The  smooth  streamline  contour  of  the  nozzles  gave  an  essentially 
uniform  exit  velocity  profile.   The  special  configuration  adopted  for 
the  control  jet  walls  prevented  the  jet  from  attaching  to  the  side  walls 
through  the  Coanda  effect.   In  addition,  two  large  holes,  as  shown  in 
Fig.  2,  were  provided  in  the  upper  and  lower  cover -plates,  in  order  to 
minimize  the  restriction  to  entrainment  to  the  resultant  jet. 

The  supply  lines  were  1=1/4  inch  diameter  Tygon  tubing,  connected 
to  the  three  inlets.   The  flow  was  straightened  via  a  one  inch  long 
honey-comb  section  before  undergoing  a  contraction  of  ten  to  one. 

Supply  pressures  were  measured  by  means  of  inclined  manometers 
through  the  use  of  wall  pressure  tabs  located  in  the  straight  part  of 
each  nozzle. 

The  top  and  bottom  plates  had  two  slots  cut  out  to  enable  the  hot- 
wire anemometer  mechanism  to  traverse  along  lines  parallel  to  the  power 
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nozzle  axis.   A  graph  paper,  graduated  in  degrees,  was  placed  under  the 
bottom  plate  for  the  purpose  of  quickly  measuring  the  jet  deflection 
angle  by  means  of  a  small  light  string  attached  to  a  needle,   A  more 
precise  measurement  of  the  jet  deflection  angle  wass  of  course ,  made 
by  means  of  the  hot-wire  anemometer. 

Air  was  supplied  to  the  system  from  a  50hp  compressor  at  approxi- 
mately 160  psi.   Air  first  passed  through  a  micro-filter  to  remove  oil 
and  dust  and  then  through  the  primary  pressure  regulator.   The  three 
jets  were  independently  supplied  from  this  regulator.   Each  jet  supply 
line  had  on  it  a  series  of  pressure  regulators,  thus  providing  indepen- 
dently regulated  flows  into  three  rotameters.   The  rotameters  were 
identical  and  each  had  a  maximum  capacity  of  17.5  scfm.  at  14.7  psia 
and  70°F.   The  actual  flow  rates  were  obtained  from 

Q  =  (rotameter  reading,  %   of  full  flow)  x  (17.5)  x  CF(P) 
where  CF(P)  represents  the  pressure  correction  factor  obtained  from  a 
chart  provided  by  the  manufacturer. 
B0    HOT-WIRE  ANEMOMETER 

A  constant  temperature  hot-wire  anemometer,  manufactured  by  the 
Thermo-Systems ,  Model  1050,  was  used  in  the  measurement  of  velocity  and 
turbulence  intensities.   The  wire  used  (wire  #1272-T1.5)  had  a  cold  re- 
sistance of  7.30  ohms  and  an  operating  resistance  of  10.95  ohms.   The 
voltage  across  the  bridge  was  measured  with  a  digital  voltmeter  accurate 
to  0.01  volts. 

A  special  probe  holder  was  constructed  and  is  shown  in  Fig.  4„   This 
holder  which  was  fitted  into  the  slots  in  the  jet  assembly,,  enabled  the 
probe  to  move  along  lines  parallel  co  the  power-nozzle  axis  in  the  mid- 
plane.   A  micrometer  mechanism  on  the  holder  provided  transverse  move- 
ment of  the  hot-wire  across  the  jet  with  an  accuracy  of  0.001  inch. 

i 
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The  hot-wire  probe  was  calibrated  with  the  Thermo-Systems  Cali- 
brator, Model  1125,  shown  in  Fig.  4.   An  inclined  manometer  was  con- 
nected to  the  reservoir  of  the  calibrator  to  measure  the  stagnation 
pressure  and  the  hot-wire  was  placed  1/16  inch  above  the  0.150  inch 
orifice  of  the  calibrator.   For  the  details  of  calibration  analysis, 
see  Appendix  VI. 
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IV„   EXPERIMENTAL  PROCEDURE 

The  experiments  were  carried  out  for  a  combination  of  three  power 
jet  velocities  and  for  four  control  jet  bias  velocities.   The  bias 
velocity  was  defined  as  (VC+VG)/2VA.   For  each  combination,  the  velo- 
city and  turbulence  distributions  were  measured  at  various  distances  from 
the  pivoting  point,  and  the  results  were  reduced  to  actual  data  through 
the  use  of  the  computerized  calibration  curve.   Most  of  the  experiments 
were  repeated  for  the  purpose  of  checking  the  reproducibility  of  the 
data.   Finally,  the  deflection  angles  were  determined  from  the  velocity 
profiles  by  connecting  the  peaks  of  each  successive  profile  with  the 
pivoting  point.   The  deflection  angles  obtained  in  this  manner  did  not 
deviate  more  than  +0.5  degrees  from  those  obtained  through  the  use  of 
the  floating  string  and  chart  paper  combination. 

The  supply  pressure  for  each  jet  was  measured,  as  previously  cited, 
by  means  of  inclined  manometers.   The  stagnation  pressure  of  the  power 
jet  was  measured  by  one  manometer,  and  the  stagnation,  as  well  as 
differential  pressure  of  the  control  jets,  through  the  use  of  two  other 
manometers. 
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V.   DISCUSSION  OF  RESULTS 

In  the  following,  first  the  experimental  results  obtained  with  the 
hot-wire  anemometer  will  be  discussed  and  then  a  comparison  will  be 
made  between  the  theoretically  and  experimentally  obtained  deflection 
angle. 

Figures  6  through  54  represent  the  velocity  and  turbulence  intensity 
distributions  as  a  function  of  x/a  and  U/VRO  in  terms  of  normalized 
control  jet  velocities,  deflection  angle,  and  the  distance  from  the 
pivoting  point.   It  is  apparent  from  all  of  these  figures  that  the  con- 
trol jets  retain  their  identity  up  to  a  distance  of  2a  from  the  pivoting 
point.   This  region  will  be  called  the  mixing  region.   For  distances 
larger  than  approximately  2a,  the  velocity  profiles  acquire  a  nearly 
Guassian  distribution.   The  turbulence  intensity  distribution  exhibits 
essentially  four  peaks:   two  in  the  outer  edges  of  the  resultant  jet 
because  of  the  mixing  of  the  jet  with  the  ambient  flow  and  two  peaks  in 
the  regions  where  the  control  jets  mix  with  the  power  jet.   As  the  dis- 
tances increase  from  the  pivoting  point,  the  turbulence  intensity  in  the 
inner  regions  decreases  and  eventually  reduces  to  a  single  minimum  at  the 
axis  of  the  jet.   It  is  therefore  apparent  that  the  application  of  the 
free-streamline  theory  to  the  mixing  of  jets  of  unequal  velocities  is 
not  entirely  valid  in  the  mixing  region  since  there  are  three  independent 
free-streamline  velocities.   In  the  fully  established  region,  the  velocity 
profile  is  Gaussian  and  thus  far  from  being  uniform.   Consequently,  the 
use  of  the  single  jet  velocity  is  not  fully  justified  unless  special  pre- 
cautions are  taken  to  make  sure  that  the  gross  continuity  of  the  mass 
flow  is  satisfied  by  properly  selecting  the  apparent  jet  width.   In  so 
doing,  the  maximum  centerline  velocity  may  be  used  as  the  free-streamline 
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velocity  and  the  effective  jet  velocity  may  be  calculated  from  the  equa- 
tion of  continuity.   The  theoretical  analysis  presented  herein  makes 
good  use  of  this  fact  and  shows  that  the  results  so  obtained  are  in  good 
agreement  with  those  obtained  experimentally. 

Figures  55  through  57  show  a  superposition  of  the  velocity  profiles 
and  the  evolution  of  the  mixing  from  three  independent  jets  to  a  single 
profile.   Similar  graphs  are  presented  in  Figs.  58  and  59  for  the 
turbulence  intensity.   Figure  60  is  a  three-dimensional  plot  of  the 
evolution  of  the  velocity  profiles  as  a  function  of  the  distance  from 

the  pivoting  point. 

The  deflection  angle  is  shown  in  Figs.  61  through  63  as  a  function 

of  the  sum  and  difference  of  the  normalized  control  jets  for  all 
Reynolds  numbers.   It  is  immediately  apparent  that  the  deflection  angle 
is  a  linear  function  of  the  normalized  differential  control  jets  velocity 
and  that  the  deflection  angle  does  not  vary  appreciably  with  the  normal- 
ized average  of  the  control  jets.  These  results  suggest  that  the 
characteristics  of  beam-deflection  amplifiers  be  evaluated  in  terms  of 
the  parameters  cited  above  rather  than  those  previously  used  by  other 
investigators.   Although  the  experimental  results  presented  herein  are 
for  equal  control  and  power  jet  nozzles,  there  Is  every  reason  to  be= 
lieve  that  similar  correlations  will  exist  between  the  jet  deflection 
angle  and  the  two  control  jets  parameters. 

A  comparison  of  the  experimental  results  with  those  obtained 
theoretically  for  a  normalized  average  control  jet  velocity  of  unity 
is  shown  in  Fig.  64.  As  previously  cited,  the  final  jet  width  d  used 
in  the  calculations  was  obtained  by  dividing  the  total  flow  rate  by  the 
maximum  jet  velocity.   It.  is  believed  that  there  is  at  present  no  other 
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way  of  applying  the  free-streamline^  theory  to  the  deflection  of  viscous 
jets  for  in  the  mixing  region  there  are  three  different  jet  velocities 
and  three  fairly  uniform. velocity  distributions  and,  in  the  fully  estab- 
lished region,  a  single  maximum  with  a  Gaussian  distribution.   Thus,  the 
only  logical  assumption  to  be  made  is  that  the  free-streamline  velocity 
used  in  the  analysis  is  equal  to  the  maximum  of  the% fully  established  pro- 
file and  that  the  effective  jet  width  is  determined  by  the  equation  of 
continuity.   In  reality,  it  appears  that  the  maximum  centerl.ine  velocity 
is  somewhat  higher  than  that  predicted  by  the  free-streamline  theory  and 
that  the  effective  jet  width  is  somewhat  smaller  than  that  used  in  the 
analysis.   Be  that  as  it  may,  no  second-order  corrections  were  applied 
to  the  analysis  in  order  to  preserve  its  generality.   The  extension  of 
the  analysis  to  other  velocity  ratios  will  be  separately  undertaken  in 
the  future. 

The  stagnation  pressures  at  supply  reservoirs  of  the  three  jets 
were  determined  as  previously  explained  and  are  plotted  in  Fig.  65. 
Mean  values  rather  than  the  actual  data  points  are  shown  in  this  graph 
because  the  data  points  fell  fairly  consistently  on  straight  lines. 
It  is  immediately  apparent  that  the  normalized  differential  control  port 
pressures  are  linear  functions  of  the  deflection  angle  and  depend  weakly 
on  the  mean  pressure  of  the  control  ports.   A  simple  analysis  of  the 
jet  deflection  through  the  use  of  the  momentum  balance  equation  shows 
that  the  deflection  angle  so  determined  is  considerably  higher  than  that 
obtained  experimentally  or  theoretically  even  when  the  pressure  differ- 
ences between  various  nozzles  are  taken  into  consideration. 
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VI  o   CONigUJSIONS 

The  foregoing  theoretical  and  experimental  study  warrants  the 
following  conclusions: 

(a)  The  jets  preserve  their  identity  in  a  region  of  approxi° 
mately  2a  from  the  pivoting  point; 

(b)  The  velocity  profiles  become  Gaussian  beyond  the  mixing 
region; 

(c)  The  free  streamline  theory  predicts  the  deflection  angles 
fairly  accurately  provided  that  an  effective  jet  width 
(calculated  in  terms  of  the  mass  flow  and  the  maximum 
centerline  velocity)  is  incorporated  into  the  analysis; 

(d)  The  previous  models  cited  in  the  introduction  do  not 
predict,,  or  are  not  capable  of  predicting;,  the  jet  de- 
flection angle  to  an  accuracy  demanded  by  the  beam  de- 
flection amplifiers;  and  finally, 

(e)  The  normalized  differential  control  reservoir  pressures 
vary  linearly  with  the  deflection  angle.   It  is  thus 
believed  that  an  analysis  such  as  the  one  presented  here, 
together  with  one  experimentally  determined  parameter,  is 
fully  capable  of  predicting  the  performance  of  proportional 
amplifiers.   The  knowledge  of  the  jet  deflection  angle  and 
the  Gaussian  shape  of  the  velocity  profile  provide  through 
a  simple  integration  the  determination  of  the  amounts  of 
flow  at  a  given  amplifier  receiver. 
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APPENDIX  I 
EVALUATION  OF  THE  INTEGRAL 
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*dv  =  t  dt 
one  has 
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The  general  solution  of  this  integral,  irrespective  of  whether  the 
radicals  are  real  or  imaginary,  is  given  in  Grobner  233.3b  [Ref.  6]  as, 
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The  constants  M  and  N  may  be  determlneds  for  the  different  regions  of 
the  t  plane,  as  follows: 
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b.      At    the   point   B, 
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c.      At   the   point        H, 
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d.   At  the  point  F, 
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Finally,  XL      reduces    to 
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Evidently,    Eq.    (28)    is   equivalent   to 
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APPENDIX  II 


EVALUATION  OF  THE  INTEGRAL 
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For  Eq.    (31), one   has 
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For  Eq.    (32)   we  have 


2=.  sin*  k. 


(35: 


ttrvv 
Qs»oo 


which  reduces  to 


V    <*>-* 


^-k1  oj  l~te1'  +Jcl-<1 


(36) 


For  Eq.    (33),   using   the   same   procedure   as   above,    one  has 

2.     i        h-$x 


Iw 


V^*-  ^  <|   \  -  VC^  '  +  «Jfl*-kV 


:37) 


For  Eq.  (34)  the  inclusion  of  the  discontinuity  at  t  s  -f  yields 


Q-^Oo 


which  after  simplification  reduces  to 


_  _*_W-£l 


^*.k-     /i^+JT^ 
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Replacing  Eqs.    (35),    (36),    (37)   and   (38)    in  Eq.    (15)    one   has 


a  +  2 


&    ^    d  n<\&    +    2  J 


\-K'- 


X 


k 


.     *fc     I, 


M       -< 


|        .    :' 


NJc'-k  *»'^     +,C-k* 


(39) 


■   +  \.  \  ■' .  <  • 
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APPENDIX  III 


EVALUATION  OF  THE  INTEGRAL 


/>+« 


*c  t- 


3 


x  w; 


Consider  the  following  typical  integral 


+\ 


'tM 


<if 


O-tf   b*p 


-i 


Multiplying  by 
have   />■+> 


f1-! 


(40) 


ce-o  db 


..^-■He-^K^p) 


(41) 


the  numerator  and  the  denominator,  we 


Dividing  (t  -1)  by  (t  +  p) 9  we  can  express  the  integral  as 
„4\  -a.1  »M 


r+* 


t  dtr 


at 


Vct'.«Htw) 


-(p-) 


dtr 


(t+pW(tu%M 


(42) 


Note  that  the  first  integral  in  Eq.  (42)  is  an  odd  function  of  t.   Thus 
its  integration  from  -1  to  +1  reduces  to  zero.   The  same  result  is 
obtained  for  the  third  integral  when  p  =  0. 

From  Grbbner  [Ref.  6]  formulas  244.8a  4  and  244.8a  7  one  has 


dtr 


and 


\|(tM(t2-e) 


dt 


•  ±  fM 


(43) 


Uk. 


"t^WCt'-oUMO"       (>^U 


(44) 
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where 

l  +  k 
and    *P     is  given  by 


2-Ct-k) 


(45) 


(46) 


and     P(f,i\)  is  defined  by 


•-«  S »'  i\  ii» 


d* 


V(i-*'Ki-KV) 


and 


*Lvf,<**h)  is  defined  as   1* 


0-otV)  -{0-ynCH^ 
In  taking  the  limit  of  Eqs.  (43)  and  (44),  it  must  be  noted  that  the  in- 
tegral from  -1  to  +1  is  to  be  performed  by  going  from  -1  to  -  90   ,  to  +  yo 
to  +1.   Replacing  the  value  of  t  in  Eq.  (46),  it  is  seen  that  ij>  goes  from 

5  to  zero  when  t  goes  from  -1  to  +1,  in  the  sense  described  above. 
Replacing  Eqs.  (43)  and  (44)  in  Eq.  (42)  and  noting  that  V   goes  from  — 
to  zero,  we  have: 


M 


t\t  dt 


Jt-e       b+ 


4k         k+p  Ck+M       L        fc+i^O      /J 


where  i\(.^J     is   the  complete  elliptic   integral  of  the   first  kind  and 

h        is   as   defined  by   Eq.    (45). 
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Replacing  Eq.  (47)  in  Eq„  (40)  and  replacing  P  by  its  appropriate .value, 


we  have 
+  »/ 


/  r 


lb 


\(u^ 


\-vck 

k-t-c 


o( 


i*-t±i 


k- 


^ 


y 


^k+f  i b 


KM 


* 


(k+»}  (.<-+»}  y 


The  singularity  at  t  =  -f  is  evaluated  separately  in  Appendix  V. 
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APPENDIX  IV 
EVALUATION  OF  COMPLETE  ELLIPTIC  INTEGRAL  OF  THE  THIRD  KIND 

In  the  general  case  of  Complete  Elliptic  Integrals  of  third  kind 


|(JU,°0)  <49> 


According  to  C.  Heuman  [Ref.  5]  they  can  be  reduced  to  combinations  of 

2 
first  and  second  kind  elliptic  integrals,  depending  on  the  value  ofc<  , 

that  is  depending  on  whether  the  integrals  are  hyperbolic  or  circular. 

In  the  present  case 


where    p   can  be  f,  c,  or  -g  and 

a.   For  the  case  of  p  ■  f,  we  have     ^  ""  >  I 

4  ia%-\}{«?-$) 


where 


f>  =.5^(ir) 

b.   For  the  case  of  p  =  c,  we  have  h  <c  ^ z <  ! 


and  the  transformation  becomes 


(50) 


(51) 


(52) 


Tf(lL   >    «*    K\    -     £        o<   AcU  >M  (53) 
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where 


=  sin  |  *  -k' 


c.   For  the  case  of  p  s  -g,  we  have     <X  <  O 
and  the  transformation  becomes t 


TT(|  ,«*  £)   .  JkX  _    TC   e<TJL(0,K)  (54) 


where 


<()  =  si  /\  .  / 2* 


*  t  t 


where  the  Heuman's  Lambda  function  —M_c   is  defined  as 


and  the  Jacobi's  Zeta  function  L.  is  defined 


as 


Z(&,^  -       ECp.h)  -  E  F(f,h)  (56) 

'  K. 

where 

Ev^/^h)  is  the  elliptic  integral  of  the  second  kind    (57) 

F  (A  ,  Vi  )  is  the  elliptic  integral  of  the  first  kind     (58) 

and        E  ~  cCa.,K)  (59) 

K-F(|A)  (60) 

and         \\      ~     I  -  \\ 

Subroutine  CEL3  was  written  to  calculate  the  complete  elliptic  integrals 

2 
of  the  third  kind  for  all  the  cases  of  <$(    "  cited  above. 
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APPENDIX  V 


EVALUATION  OF  THE  DISCONTINUITY  OF  THE  INTEGRAL 


6fc 


(61) 


MWtt^)(e-.) 

Since  there  is  a  discontinuity  at  t  =  -f,  we  have  to  go  to  the  limiting 
process  as  follows 


r< 


f+£  +» 


\^(t)     -_        1  fCt)     +        |W  <"    £*°  (62) 


As   calculated   in  Appendix   IV 


dfc 


'MWo^GrV)      c 


_^L_ .  i  KK)  -  ^Tffa*^^     (63) 


where 

siiff  c    ^-<X^)  (64) 

When  taking    f(^)   dr    »  we  have  J  t(f)     which  means  that  *~(".  '0 

reduces  to  K  (complete  elliptic  integral  of  the  first  kind)  and  TT  (  |  pVL^\\  \ 

2 

reduces  to    fl  (  ^  ,  ||  )   (complete  elliptic  integral  of  the  third  kind). 

The  only  term  which  has  a  discontinuity  is  '{O^i'i  )  since  !~(f  i\\  is  continuous, 

At  ta-jt£. 

and   at       t  *-£-£. 

t_*     sin"     lO-HQl'+KO 

Then  Eq.  (62)  becomes 

Bto)-'!XH^  -v  BOM-BCD  (65) 
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Let  us  examine  in  particulars,  "BOn  )  ~"&(tz.)  t   or 

T rans forming ^[pP f <*  ,  hj  into  a  combination  of  first  kind  elliptict 
integrals  and  Zeta  and  Theta  functions9  according  to  [Kef„  4],  we 


f66' 


(67) 


where 


«^4  =    J-  iv  Sinl^C^  +  TJ  \         m    sm2m'Vs»rl,2majCM 


sin[coC«v)-  v] 


rn^\ 


1^    5vrtV\  2xr\ 


? 


2K 


-2p 


V-  e 

when  we    take      *&W\)  ~£>ll2.)     where    BC^J  is  given  by  Eq.    (67),   we   are 
left  with  B^-U^=-     MMk^lf   ^^)1  (68) 

Examining  it^,      further  we    find   that  every   term  is   independent   of  £.    ex- 
cept   Hr  ,   which   is  a   function  of    i      .      Thus  Eq»    (68)  may  be  written  as 


BU>  )-■&(%)  -X 


^(i4f((i'h^fW"^] 
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Simplifying,  Eq.  (69)  reduces  to 


in  the  limit  as  ^O,  %   "*il  "~*  X 

The  ratio  inside  the  U\,  will  reduce  to  unity  unless  F(A,  \\   )   =  +  F(  vj-  ,^) 

Replacing  A  by 


and  If   by 


Then         A-  ^ 


»?  »    iiV? 


1 


2-O-H) 


Thus  we  have  an  undetermined  case  of  the  type   A*v 


L    o 


Defining 


V2. 


=  b 


/Sk^i 


2.K 


and  using  L'Hopital's  rule,  we  have 


(70) 


which  is  equal   to 

,.       lm     b     ivy,    H>jUf    -*0 

sm  (a  £  -   ft  £,  J 


L-»o 
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Using  L'Hopital's  rule 


or 


i- 


From  Hildebrand  (Adv.  Calc.  for  Eng.)  page  352,  one  has 


where 


£t^*    |C£)  = 


Then 


wop,)   =  ,r^o     Lkii 


^L. 


■ 


and 


^  si*  If  J 

~  —  > 

J¥ 

Now  replacing  Eqs» 

(72), 

(73), 

£.  -*•    0S    one  has 

*-      iCf^^-^0  z/fu~+I 


(\-*0 


-tH*-0*Vf^Z 


and    (75)    in  Eq, 


(71) 


(72) 


(73) 


(74) 


75: 


and   letting 
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APPENDIX  VI 

DERIVATION  OF  TURBULENCE  -  INTENSITY  EQUATION 

On  a  hot-wire  anemometer  the  heat  transfer  equilibrium  equation  is 
given  by 

^  **     .    /W  B  VTT 

*w  -  "£cj  (76) 

where     I  =  current  through  the  wire 

Rw  =  operating  resistance 

K^  =  cold  resistance 

A  And  B  =*  constants 

U  =  instantaneous  velocity  of  the  fluid 

From  (76)  at  (J=   0   I  =  I0 
So 


2  2      2 

Multiplying  both  sides  of  the  equation  by  R  and  replacing  I  R  by  t   we 

w 

obtain  »  , — 


or 


e*--  E0r  -v  Cvu 


(77) 


Since  the  wire  is  kept  at  a  constant  temperature  R  and  R  are  constants. 

w      g 

Taking  the  derivative  with  respect  to  time,  one  has 

2edE  a  c£y    du  (78) 

dt    2  U   c*c 
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From  Eq.  (77)        C  >/ U   -  E"*-£0 


Then  Eq.  (78)  becomes 


where     dU  -   i/^£    ~  velocity  variation 
and       43  z.    y  Q  ~   voltage  variation 


Finally  we  have 

in  other  words,  turbulence  intensity  is  equal  to 
times  RMS  reading. 


4E 


R**  E^ 
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TABLE   I 


Corresponding  Value 

Corresponding  Value 

Point  in  the 

in  the 

in  the 

z-plane 

t-plane 

-plane 

A 

0 

V,     2. 
\   C 

B 

-k 

C 

-c 

D 

-1 

c  +  to 

E 

f 

0  4  t  <6 

F 

1 

0  +  i.  TC 

G 

8 

■  J* 

H 

i 


k 
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Figure    la     z-plane 
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Figure  1  Transformations  Planes 
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Figure    2     Jet  Assembly 
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Figure   3     Jet  Assembly 
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Figure  5   General  Setup  of  Experiment 
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FIGURE   6        VELOCITY    PROFILE  AND   TURBULENCE    INTENSITY 
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FIGURE   7      VELOCITY   PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE   8      VELOCITY   PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE  9   VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE  11   VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE    13      VELOCITY    PROFILE  AND   TURBULENCE    INTENSITY 
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FIGURE  14   VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE    18      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE    19      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE    21      VELOCITY   PROFILE  AND  TURBULENCE   INTENSITY 
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FIGURE    22      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE   25      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE   26      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE   27      VELOCITY   PROFILE  AND  TURBULENCE   INTENSITY 
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FIGURE    29      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE   30      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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5  FIGURE  31   VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE   32      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE  33   VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE   34      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE      35      VELOCITY    PROFILE  AND  TURBULENCE   INTENSITY 
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FIGURE  36  VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE   37      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE   38      VELOCITY   PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE   39      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE   40     VELOCITY    PROFILE  AND  TURBULENCE   INTENSITY 
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FIGURE  41      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE  42   VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 


93 


FIGURE  43   VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE  44  VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE  45      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE  46  VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE  48   VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE  49  VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE   50      VELOCITY    PROFILE  AND  TURBULENCE    INTENSITY 
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FIGURE  52  VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 
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FIGURE   53      VELOCITY   PROFILE  AND  TURBULENCE   INTENSITY 
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FIGURE  55  VELOCITY  PROFILES 
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FIGURE   56      VELOCITY    PROFILES 
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FIGURE   57      VELOCITY   PROFILES 
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FIGURE   58      TURBULENCE   INTENSITY 
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FIGURE   59      TURBULENCE   INTENSITY 
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FIGURE  60  VELOCITY  PROFILES  AT  DIFFERENT  STATIONS  DOWNSTREAM 
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